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We study transport in topologically disordered networks that are subjected to an environment that 
induces classical diffusion. The dynamics is phenomenologically described within the framework of 
the recently introduced quantum stochastic walk, allowing to study the crossover between coherent 
transport and purely classical diffusion. We find that the coupling to the environment removes all 
effects of localisation and quickly leads to classical transport. Furthermore, we find that on the level 
of the transport efficiency, the system can be well described by reducing it to a two-node network 
(a dimer). 
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I. INTRODUCTION 

Recent experimental advances in ultra-cold Rydberg 
gases allow for a precise control of its constituent atoms, 
including the strength of their interactions by the use of 
Forster resonances At these temperatures the ther- 
mal energy is much smaller than the interaction energies 
and the interaction dynamics is much faster than the 
displacement due to thermal motion. The dynamics is 
therefore completely driven by Rydberg-Rydberg inter- 
actions [2]. In experiments, one usually excites one atom 
into a Rydberg state that is resonantly coupled to a lower 
state. Resonant energy transfer now leads to quick coher- 
ent hopping of the excitation to nearby atoms Q • Since 
the large dipole moment of this Rydberg state leads to 
the dipole-blockade mechanism preventing the appear- 
ance of multiple Rydberg excitations in a nearby area 
[H , one can study excitonic transport in this system by 
considering a network of coupled two-level systems 

Another important example of resonant energy trans- 
fer is the study of exciton transport @ in light harvesting 
systems found in, for example, marine algae 01. Af- 
ter the absorption of a photon, the exciton is transfered 
along a series of bacterio-chlorophylls (BChl) to the reac- 
tion center. Recent experiments indicate that this trans- 
port shows coherent features, even at room temperatures 

For both of these systems it is important to study 
the effects of environmental interactions on the coher- 
ent transport of excitations [13[ . These interactions will 
in general lead to decoherence and are therefore usually 
detrimental to the transport efficiency. There are situa- 
tions however, where the presence of onsite energetic dis- 
order and decoherence rates that are proportional to the 
intra-site couplings may actually lead to faster transport 
[lij . This is for example demonstrated in [l5| for a fully 
connected network and in [r| [l7j for the FMO-complex 
found in light harvesting systems. In these works the 
environment is assumed to be weakly coupled to the sys- 
tem, leading to Markovian transport. 

On the other hand, one can also take the following 



complementary viewpoint: instead of modelling the ex- 
act type of environmental interactions by phenomenolog- 
ical master equations, one can also design a system-bath 
coupling to engineer systems that exhibit efficient trans- 
port. This idea was for example suggested in [l8| and 
[19j where such a procedure for the use of quantum state 
engineering has been proposed. 

Motivated by these viewpoints we consider systems for 
which the environmental interactions induce incoherent 
transfer of the populations. These interactions have the 
advantage that for large couplings to the environment the 
system does not run into the quantum Zeno effect, but in- 
stead behaves as if it is governed by classical diffusion. By 
suitably engineering the system and the coupling to the 
environment, we can therefore study the transition from 
purely coherent dynamics, which is usually described by 
the continuous-time quantum walk (2^, [21j , to purely in- 
coherent dynamics that is described by classical diffusion 
(or the continuous-time random walk) [2^ . 

For the mathematical framework underlying this tran- 
sition,^ use the recently introduced quantum stochastic 
walk [23| . This is a generalisation of the continuous-time 
quantum walk by also allowing for incoherent transfer 
between the sites of the system. This model also shows 
similarities to early approaches like the Haken-Strobl- 
Reineker model [HI HE]- To study the transport efficiency 
we furthermore connect our system with a source and a 
drain and model the transitions from and to the system 
by irreversible incoherent transfer. This idea was already 
implemented by several groups [TH, [H, [26l - [28j | . 

In this paper we take our model system to be a topolog- 
ically disordered network with long-range dipole-dipole 
interactions, resembling for instance a gas of ultra-cold 
Rydberg atoms. We study the crossover between quan- 
tum and classical transport and the transport behaviour 
as a function of the source and drain strengths. We find 
that on the level of the transport efficiency, we can effec- 
tively reduce our network to a dimer. 
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II. THE QUANTUM STOCHASTIC WALK 

In this section we review the mathematical details of 
both the continuous-time random walk (CTRW) and the 
continuous-time quantum walk (CTQW) and we intro- 
duce the concept of the quantum stochastic walk (QSW) 
that describes the transition between the CTRW and the 
CTQW. 



which describe the overlap of the initial state \l) with 
the final state fc) after a time t. Equivalently, one can 
also formulate the CTQW by specifying the evolution 
of the density operator p(t), which is described by the 
Liouville-von Neumann equation: 



dp(t) 
dt 



[H,p(t)] 



(5) 



Random walk 



3. Quantum stochastic walk 



Consider a network consisting of N nodes. One can 
specify the connections in the network by the connectiv- 
ity matrix A. This N x N matrix is defined by: 



fj for k = j 
Akj = \ —1 if fc and j are connected 
else 



(1) 



where fj is the functionality of the node j, i.e. the num- 
ber of bonds emanating from j. To each node fc a vector 
|fc) is associated such that the collection of all these vec- 
tors forms an orthonormal basis of TV-dimensional vector 
space. 

The dynamics on the network is governed by the trans- 
fer matrix T, which is the matrix of transition rates per 
unit time. For the simplest case, where all the transition 
rates are equal to, say 7, T is related to the connectivity 
matrix A by T = — 7 A. After assuming that the trans- 
port is described by a Markovian process, one arrives at 
the master equation [22] 



dpki (t) 
dt 



N 



i=i 



(2) 



where p k i (t) is the probability of being at node fc after 
a time t, under the constraint that one starts at t = 
in some node I. This master equation is the defining 
evolution equation for a CTRW. 



2. Quantum walk 



Now, if one places the system in an external envi- 
ronment, the full Hamiltonian takes the form H to t = 
H + H^; + Hi n t , where H is the Hamiltonian of the net- 
work, is the Hamiltonian of the environment and 
Hi n t specifies the interactions between the network and 
the environment. When the environmental correlation 
time is small compared to the relaxation time of the sys- 
tem, one can employ the Born-Markov approximation. 
This approximation results in the following general form 
for the evolution equation of the reduced density operator 
pit) of the system [29l . |30| , which is called the Lindblad 
equation: 



dp(t) 
dt 



N 



= -1 [H, p{t)\ + *kiV(L kl ,p(t)), (6) 



fc,;=i 



with the constants X k i > for all fc and / and 
V{Uup{t)) = L w p(t)L+ w - \ {hlL kV p(t)} 



(7) 



The operators L/./ are called Lindblad operators and they 
form an orthonormal basis for the space of operators act- 
ing on the system's Hilbert space. The constants Xki 
are related to certain correlation functions of the envi- 
ronment and they play the role of relaxation rates for 
different quantum channels. For notational clarity we 
introduce the dissipator D(p(t)) as 



N 



V(p(t))= X kl V(L kh p(t)). 



(8) 



fej=i 



The CTQW can be formulated in a similar fashion. 
One can interpret the basis vectors |fc) as forming a ba- 
sis for the whole accessible Hilbert space. The main 
idea is now to identify the CTQW Hamiltonian H with 
the CTRW transfer matrix T by setting H = T [20|. 
The dynamics of the system is then described by the 
Schrodinger equation 



— a u {t) = -i^H kj aji{t), 
3 

for the transition amplitudes 

a kl (t) = (k\e- mt \l) , 



(3) 



(4) 



When it is clear from the context, the explicit dependence 
on p(t) will be omitted. 

The Lindblad equation is the defining equation for the 
quantum stochastic walk (QSW) [23[. It contains both 
the coherent evolution due to the Hamiltonian and the 
incoherent evolution due to the environmental interac- 
tions. Until now we have not made any assumptions on 
the particular details of the coupling to the environment. 
Here we focus on the transitions between the CTQW and 
the CTRW, thus we consider environmental interactions 
that will eventually lead to classical diffusion. This can 
be achieved by a proper choice of the Lindblad operators 
h k i , as is shown in appendix [X] The dissipator corre- 
sponding to these Lindblad operators describes a CTRW 
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for the populations (the diagonal elements of pit)) and 
pure decoherence leading to exponentially decaying co- 
herences (the off diagonal elements of pit)). Furthermore, 
we choose the coupling constants Xki to be equal to the 
absolute value of the corresponding matrix elements of 
the Hamiltonian, that is Xki = \Hki\ — \Tu\- This en- 
sures that the CTRW has the same transfer rates as the 
CTQW. 

The full evolution equation of our system is then given 
by a linear combination of the CTQW and the CTRW 
combined with pure decoherence, together with a scaling 
parameter a |23j : 



dpjt) 
dt 



(1 - a)£ coh (p(t)) + aVipit)), a e [0, 1], (9) 



with £ C oh(p(^)) = — i [Hs,p(i)] representing the genera- 
tor for purely coherent transport. In the limit a — > we 
obtain the CTQW and in the limit «->lwe obtain the 
CTRW. 



III. SOURCES AND DRAINS 
A. General considerations 

In order to incorporate sources and drains into the 
system, we consider the following scenario: for a sin- 
gle source we augment the network of N nodes by an 
additional node |0) = (source). In order to prevent the 
excitation from flowing back into the source this node will 
be incoherently coupled to S nodes of the system, with 
1 < S < N. The incoherent nature of the coupling im- 
plies that we do not couple the source to the network by 
the Hamiltonian, but that we use the Lindblad formalism 
to describe an incoherent hopping from the source to the 
network, i.e. = for all the nodes |fc) comprising the 
original network. 

Similarly, we include a drain as an extra node 
\N + 1) = | drain) that is incoherently coupled to D nodes 
of the system, with 1 < D < N. Thus the total dimen- 
sion of the systems Hilbert space will be N + 2 and the 
reduced density matrix of the source-network-drain sys- 
tem will be a (N + 2) x (N + 2) matrix. 

When the source (drain) is coupled to more than one 
node of the network, there arc different ways to model the 
transition from the source (drain) to these nodes and the 
resulting dynamics can vary strongly. However, in this 
paper we only consider the simplest case where both the 
source and drain are connected to only one node, leaving 
the more general case for future work [3lj 

The source node |0) will be coupled to node \k) of the 
network by introducing the Lindblad operator likfi — 
\k) (0| (see e.g. [Ill H3]). We will denote its coupling 
strength by IV Similarly, the drain node \N + 1) will be 
coupled to node \l) by introducing the Lindblad operator 
Ljv+U = \N + 1) {l\ and its respective coupling strength 
will be denoted by 7;. 



We pause to note that there is a difference between the 
Ffc and ji introduced here and in previous work |32L [33[ . 
There a drain connected to a node m) was modelled 
by considering the effective Hamiltonian H e // = H — 
*7m \fn) (m\. As is shown in appendix [Bl in the Lindblad 
approach this is equivalent to a coupling strength equal 
to 27,„. 

The dynamics of the reduced density operator pit) 
is then given by the usual Lindblad equation for the 
quantum stochastic walk plus the extra dissipators corre- 
sponding to the Lindblad operators that model the cou- 
pling of the source and drain to the network: 



dpjt) 
dt 
with 



(l-a)£ coh (p(t))+aD( / o(t))+r s+d (p(t)), (10) 



C s +d(p(t)) =r fc X>(L fc ,o,p(*)) + 

rf,V(L N+1 ,,,p(t)). (11) 

Note that both £ co h and £>(p(i)) only act on the subspace 
spanned by the network nodes, while acts on the 

combined source-network-drain system. 

As an initial condition we will always choose to start 
in the source node, so p(0) = |0) (0|. By the explicit 
form of the dissipators given by Eq. (|A1|) and by using 
that Hok — HN+i : k = for k = 0, . . . N + 1, we have 
the following expressions for the coherences between the 
source resp. drain and the rest of the network: 



dpjojt) 
dt 

dpj,N+iit) 



= -^r fc p i0 (t), vj^o 
1 



(12) 



= -- llPh N+1 (t), + (13) 

Our choice of initial conditions then implies that all these 
coherences are identically zero. This means that the den- 
sity operator pit) can be written in the following block- 
diagonal matrix form: 



Pit) 



Pooit) 
pit) 







PN+l,N+lit) 



(14) 



with pit) the corresponding density matrix restricted to 
the subspace spanned by the nodes of the original net- 
work. 

Although we cannot make any detailed predictions for 
the dynamics of our system, we can make a general state- 
ment on the temporal behavior of the population of the 
source node. From Eqs. (flU)) and (fTTj) it follows that 



leading to 



dpoojt) 
dt 



Pooit) 



-T k pooit), 



(15) 



(16) 



Thus, the population of the source will in any case decay 
exponentially. Since the source is incoherently coupled to 
the network, eventually all its population will be trans- 
ferred from the source to the network. 
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FIG. 1. Schematic representation of a single realisation of a 
topologically disordered network in a sphere with long-range 
dipole-dipole interactions. 



B. Transport efficiency 

In order to gain insight into the efficiency of the trans- 
port from the source to the drain, we adopt the following 
definition of the transport efficiency which is similar to 
the one used by Cao and Silbey [T3] and by Aspuru-Guzik 
et. al. 17]: 77(a) is the expected survival time (EST), i.e. 
the expected amount of time that the excitation will re- 
main inside the source and the network nodes: 

N 



](a) = / dt^2p kk (t,a) 
Ja k=o 



dt[l - PN+i,N+i(t, a)} . (17) 



We stop to remark that the contribution of the source 
node to Eq. (IT7|) amounts to a constant, since the decay 
of poo(t) is exponential, see Eq. (ITtH) . If 77(a) is a small 
number this means that transport to the drain is rela- 
tively fast, and vice-versa. 77(a) can be calculated from 
the Laplace transforms of the populations p k k (t) , see ap- 
pendix. [C] For instance, for a network consisting of only 
one node |1) that is connected to a source and an drain, 
the EST 77(a) can be easily computed, see appendix [El 



77(a) = i/r + 1/7, 



(18) 



where T and 7 are the source resp. drain strengths. The 
EST is therefore independent of a and the transport be- 
comes more efficient with increasing values of 7 and T. 



to the sphere. Previously, a similar system was studied, 
but without any coupling to the environment [34]. Also 
in [35|, |36[ a similar set-up was considered, but there the 
focus was not on the quantum-to-classical crossover we 
study in this paper. 



A. The model 

We consider a network of N sites located inside a 
sphere of radius R. At both sides of the sphere we place 
a node and they are denoted by |1) resp. \N). The other 
N — 2 nodes are chosen randomly. Furthermore, we con- 
nect a source |0) to |1) and a drain \N + 1) to \N); their 
respective coupling strengths are given by T and 7, see 
Fig. [1] for an illustration of this system. As was noted 
before, we choose dipole-dipole interactions between the 
nodes, decaying as d^ 3 with dki being the distance be- 
tween the nodes fc and I. The matrix elements of the 
Hamiltonian then take the following form: 



(*|H|/) 



-d~ 3 

u kl 



iovk^l 



(19) 



In this particular case, the QSW equation descrbing the 
transition from the CTQW is given by Eq. (flU]) . 



dp(t) 
dt 



N 



= (1 — a)£ co h + a ]T \U kl \V(L kl ) + C s+d , (20) 



with the source-drain superoperator 

C s+d (p(t))=TV(L 1 , ,p(t))+ 1 V(L N+hNl p(t)). (21) 

For each realisation r of the system, with r = 1, ... ,72., 
we can numerically solve this equation and find the solu- 
tion p( r ) (t) . We calculate the ensemble average 



p(*) = ^Ep (r) w 



(22) 



in order to obtain a global picture of the general dynamics 
of this system. 



B. Numerical results 



IV. TOPOLOGICALLY DISORDERED 
NETWORKS 

We now apply the above concepts to systems that ex- 
hibit topological disorder and long-range dipole-dipole in- 
teractions. Here we model these systems by considering a 
random configuration of nodes in some bounded region of 
space, say a sphere, which is coupled to an environment 
that induces diffusive behaviour. To study the transport 
efficiency we furthermore connect a source and a drain 



For all numerical results that will be discussed in 
this section, we assume that N = 7 for computational 
reasons. For the total number of realisations we take 
1Z — 4000. Without loss of generality we can take for the 
sphere unit radius, since other radii would merely rescale 
the average distances between the nodes and therefore 
would only provide a rescaling of the time axis. The 
precise value of the source strength T is not particulary 
important since it only effects the flow rate into the net- 
work, but does not qualitatively influence the flow to the 
drain: its contribution to the EST 77(a) is just a term 
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(a) a = (b) a = 




FIG. 2. Populations of a topologically disordered network with N — 7 nodes for a = 0, (a) and (b), and a = 0.01, (c) and 
(d), with r = 0.5 and 7 = 1. Figures (a) and (c) show the populations of node |1) and node |iV) as well as the sum of the 
populations of the remaining N — 2 nodes of the network. Figures (b) and (d) show the populations of the source and of the 
drain. 



equal to 1/r. For the rest of this section we therefore 
choose T = 0.5. 

As was shown in Eq. (|16l) . the initial excitation will 
decay exponentially with rate T from the source node into 
the connected node of the sphere, after which it flows to 
the drain site. In Fig. [5] we show the populations of 
the system for the two different values of a = and 
a = 0.01. Here we indeed observe an exponential decay 
of the population of the source (see poo(t) m Figs. Hfb) 
and[2]Jd)), and an initially high population of node |1) 
(see pu(t) in Figs. [5Ja) andrjfc)). 

An interesting observation arises from this figure: in 
the case of purely quantum mechanical transport (a = 0) 
we observe a form of short-to-intermediate time localiza- 
tion, with about 23% of the population remaining in the 
network after t = 40, see Figs. HJa), Hfb) and Fig. 
. Usually however, as time progresses there is still some 
population leaking into the drain. In the ensemble aver- 
age eventually everything is then transferred to the drain. 
After switching on the environmental induced diffusion 



the localization effect completely vanishes, as can be seen 
in Figs. He) andHd). 

In Fig. [3] we show the total population of the source- 
network system (excluding the drain) for different values 
of a (note the log-lin scale). For a — we observe, 
in the inset of Fig. [3] that shows the curve in a log- 
log scale, a power-law decay for intermediate times with 
77(a) ~ t~P and (3 ~ 0.21, which is characteristic for 
quantum walks in topologically disordered networks [34| . 
Already for small values of a, however, this behaviour 
quickly vanishes, ultimately leading to a pure exponential 
decay for a = 1. Here we find that r]{a) ~ e~ M * with 
fi = 0.247. This means that classical diffusion already 
dominates over quantum transport for relatively small 
values of the environmental interaction. This conclusion 
holds true also for other values of N (calculations which 
we do not show here). 

Figure H] shows a plot of the EST 77(a) for this system 
for various values of 7. As was already noted before, 
the transport efficiency is positively influenced by the 
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FIG. 3. Total population of the source-network nodes for 
various values of a = 0, 10" 4 , 10" 3 , 10" 2 , 10" 1 , 1, V = 0.5 and 
7 = 1. The inset displays the curve for a = in a log-log 
scale for the region that shows powerlaw decay. 
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FIG. 4. The EST 77(a) as a function of a for T = 0.5 and for 
different values of 7. 



environmental diffusion, since this removes the effects of 
localisation. We see that this is indeed also reflected in 
the computed transport efficiency, namely all curves for 
77(a) show a monotonous decay with increasing a. We 
also observe that higher values of 7 lead to similar curves 
for the EST as a function of a, but with a lower overall 
value. Therefore larger trapping rates and larger values 
of a lead to faster transport to the drain. 

This leads us to conjecture that for systems with 
quenched disorder and long-range (dipole-dipole) inter- 
actions transport will be on average dominated by clas- 
sical diffusion, provided that the system is coupled to an 
environment that induces diffusive behavior. 



V. DESCRIBING THE DISORDERED 
NETWORK WITH A DIMER 

In this section we investigate to what extent the topo- 
logically disordered network of the previous section can 
be described by an effective dimer model. To do this, 
we first focus on the properties of the EST 77(a) for the 
general case of a heterodimer that is coupled to a source 
and drain. Its Hamiltonian is defined by: 



H 



-V 
-V A 



(23) 



where V represents the hopping rate between the nodes 
and A represents the energetic disorder between the two 
nodes. 

In contrast to the case of the monomer (see Appendix 
|E|) , there exist multiple ways of connecting the source and 
drain to the dimer. Here we focus only on the following 
case: we connect the source to node |1) and the drain to 
node 1 2). Only this configuration allows us to describe 
the topologically disordered network with a dimer. 

The master equation (|10p can be solved analytically in 
this case. However, the full expression does not provide 
much insight, but leads to the exact expression for the 
EST 77(a), see Appendix |Dj In the limit A — > we 
obtain: 



77(a) 



2 

1 



1 

r 



1 

Va 



4(1- a) 2 



aV(4 - 8a + 6a 2 ) + a 2 7 ' 
For 7 = r we find the following two limiting cases: 

7 2 + 4A 2 + 3 



lim 77(a) 

a— 



4I/ 2 7 



, X 3 1 

km 77(a) = - + — 

a— >1 7 V 



(24) 

(25) 
(26) 



For purely classical diffusion the transport efficiency 
is therefore completely determined by the hopping rate 
V and the source and drain rates 7. This is however 
not the case for quantum mechanical transport where A 
has a large influence on the transport efficiency. Larger 
values of A lead to a quick increase in 77(a), as can be 
observed from Fig. [B]and Eq. (1231) . 

In Fig. [5] we show the EST 77(a) for the case A = 

and for various values of the drain strength 7. We 
observe that higher values of 7 lead to lower values of 
the EST and therefore to faster transport. But unlike 
the disordered network, we don't observe a monotonous 
decay with increasing a, but instead a maximum around 
a = 0.77. 

The more general case of A ^ is shown as a contour 
plot in Fig. [5] Here we observe that the EST is positive 
for all values of A and that it approaches the constant 
value 1/V + 3/7 in the limit a — > 1. For values of |A| > 

1 we do not observe a maximum of the EST anymore 
but instead a monotonous decay with increasing values 
of a, resembling the EST for the topologically disordered 
network, see Fig. 2] 
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FIG. 5. The EST 77(a) of a dimer for various values of 7, with 
V = 1, A = and V = 0.5. 
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FIG. 6. Contour plot of the EST 77(a) for a dimer as a func- 
tion of a and the energy disorder A for 7 = V — 0.5 and 
V — 1. The distance between the contour levels is 0.5 and the 
more outward lying contours have a higher value of 77(a). We 
observe that 77(a) is also positive for all values of A. Further- 
more, the EST approaches the constant value 1/V + 3/7 = 7 
in the limit a — ► 1. 



After having studied the general properties of the EST 
of a dimer, we now proceed with fitting it to the EST 
of the topologically disordered network. In Fig. [7] we 
provide a fit of the EST for the disordered network with 
r = 7 = 0.5. The corresponding parameters for the 
dimer are 7 d = 1.23, T d = 0.19, V = 0.61 and A = 1.8, 
where the subscripts d refer to the source and drain 
strengths of the dimer. In performing this fit, we have 
only fitted T^, 7^ and V because we have the freedom 
to choose the disorder strength A, as long as |A| > 1. 
Higher values of A then correspond to larger values of 
both jd, Td and V. This can be easily understood since 
in order to overcome a higher energy barrier one must 



FIG. 7. Fit of the EST of the topologically disordered network 
with 7 = r = 0.5 to that of a dimer. The fitting parameters 
are given by: -y d = 1.23, r d = 0.19, V = 0.61 and A = 
1.8. Similar fits can be made for other values of 7 and V. 
One observes that the transport efficiency of a topologically 
disordered network can be well described with a dimer. 

increase the source strength and at the same time also in- 
crease the drain strength. For other values of 7 and T for 
the disordered network, the fit behaves in a similar fash- 
ion. Therefore, the transport efficiency of a topologically 
disordered network can be well described by modelling 
the system as a dimer with finite energy disorder. 



VI. SUMMARY 

In this paper we studied the transport efficiency of an 
excitation moving from a source via a network to a drain. 
As a model of many interesting physical systems such as 
ultra-cold Rydberg gases, we chose a topologically dis- 
ordered network with long-range interactions of dipole- 
dipole type. In particular, we focused on the crossover 
between purely quantum mechanical transport and en- 
vironmentally induced diffusion which we phenomeno- 
logically modeled by employing the recently introduced 
quantum stochastic walk. 

Without any environmentally induced diffusion we ob- 
served localization at short to intermediate times of the 
excitation. After switching on the environment this ef- 
fect quickly vanished, with the total population in the 
network decaying exponentially to the drain. This ef- 
fect is furthermore largely indepedent on the number of 
nodes in the network. We can therefore conclude that in 
a system with randomly placed nodes and dipole-dipole 
interactions between the nodes, transport is on average 
mainly dominated by classical diffusion, provided that 
the system is coupled to an environment that induces 
diffusive transport. 

As a measure for the transport efficiency we used the 
expected survival time (EST) 77(a), defined to be the 
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expected amount of time that it takes to transfer the 
excitation from the source to the drain. We found that 
for all values of the source and drain strengths, transport 
is most efficient in the purely classical case. Furthermore, 
we showed that on the level of the EST the topologically 
disordered system can be mapped on a dimer with a finite 
energy difference between the nodes. 
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Appendix A: Random walks in the Lindblad 
formalism 



Here we provide a short derivation on how to obtain a 
random walk in the Lindblad formalism. In order to pro- 
vide this connection, we consider Lindblad operators of 
the form h k i = \k) (l\ for k, I = 1, . . . N. Since the Lind- 
blad operators form a complete orthonormal basis, we 
can expand the dissipators T>(Lki, p) in these operators: 



1 



£>(Lm) = piil>kk - Pii^u — 2 Y (pij^ij + Pji^n) ■ 

(Al) 

Suppose now that we reparametrize the Lindblad equa- 
tion as in Eq. © and take a = 1. The dynamics is then 
completely governed by the dissipator T>(p). Inserting 
the above expansion in the Lindblad equation results in 
the following differential equations for the diagonal com- 
ponents of the density matrix: 



a N 

~T, — 2_^\*klPll — MkPkk), 
' 2 = 1 



and similarly for the off diagonal components: 

dpkm 



dt 



1 - 

-Y 

2 ^ 



(Ajfe + Xjm)Pk 



(A2) 



(A3) 



This implies that the equations for the coherences pkm 
decouple from the populations pkk ■ Furthermore, we can 
rewrite Eq. (|A2[) for the populations in the form of Eq. 
©, with 

N 

T k i = (1 - 5 k ,i) Xki + S k ,i Y X i- k - ( A4 ) 

The matrix T formed by the coefficients T k \ can then 
be interpreted as the transfer matrix for a CTRW. The 
coupling constants X k i therefore represent the transfer 
rates of moving from node \l) to node 



Appendix B: Relation to modelling traps with 
effective Hamiltonians 



We now provide a short proof of the equivalence of 
modelling a drain with Lindblad operators and of mod- 
elling a drain by introducing an imaginary part to the 
Hamiltonian, as was done in previous works [32, 33] . Sup- 
pose that we connect a drain to node \N) with strength 
7n. One can then introduce the effective Hamiltonian 
H e // = H - if with f 1 = 7jv \N) (N\. This leads to the 
following modification to the Liouville equation: 



dps . r„ 1 

- w = -t[ll,ps] 



(Bl) 



Expanding the last term of this equation in terms of the 
Lindblad operators results in 



1 N ^ 

Pnn^nn + 2 Y^ (pNj^Nj 

3 = 1 



(B2) 

These are clearly the same terms as the one that would 
be obtained from Eqns. (IA1|) in combination with Eq. 
(fTTjl . after projecting the dissipator onto the subspace 
spanned by the network nodes. The only difference is 
the factor of two that arises in the coupling strength 7jv, 
but this is a matter of convention. 



Appendix C: Computation of 77(a) with a Laplace 
transform of the density matrix 

The EST 77(a) can also be written in terms of the 
Laplace transforms p kk (s) of the diagonal components 
of the density matrix: 



77(a) = limV/ dte~ st Pkk (t) (CI) 

N 

= hmV/5 fefe ( S ). (C2) 

S— 5-0 * * 



fe=0 



Writing the components of the density matrix in a vector 
p(t) = (pa (t),pn(t), . . .,p N+ltN+ i(t)) leads to 



dp{t) 
dt 



= Cp(t)> 



(C3) 



with C being the superoperator including all coherent, 
diffusive and incoherent terms. This then results in the 
following equation for their Laplace transforms: 



(sI-C)p(s)=p(0)- 



(C4) 
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Appendix D: Complete expressions for the EST of 
the dimer 



For the dimer with a nonzero value of the energy offset 
A we find the following expression for EST 77(a): 



17(a) 



1 



1 

a 



1 

V 



/(«) 



Eqns. (fTU)) - (fl~3")) , we obtain a master equation where all 
non-diagonal elements are decoupled, each of which ac- 
quires a simple exponential solution. Our choice of the 
initial condition p(0) = |0) (0| thus yields Pkj{t) = for 
k 7^ j. For the diagonal elements we find the following 
equations: 



3(a) 



(Dl) 



7 r 

The function /(a) is given by: 

/(a) = 4(1 - a 2 ){2Va + 7 + aA), 
and the function g(a) is given by: 

g(a) = W 2 a(2 + a(3a - 4)) + 4V(1 + 2(a - l)a)(7 + aA) 



(D2) 



+a(7 2 + 2a7A + (4 + a(5a - 8))A 2 ) 



(D3) 



Note that both /(a) and g(a) are independent of the 
source strength T. 



Appendix E: A monomer with source and drain 

Here we provide a simple example to illustrate the con- 
cepts introduced in section Hill We take the network to 
be only a single node |1) such that the source-network- 
drain system consists of three nodes in total. For nota- 
tional simplicity we will write F = Ti and 7 = 71. The 
Hamiltonian H is proportional to |1) According to 



poo(t) = -I>oo(<) 

pn(t) = rpoo(*) - ipn(t) 

P22(t) = 7Pll(*)- 

which have the solutions 

Poo(t) = exp(-rt) 

Pn(f) = fT~ — (exp (-7*) - exp (-Ft)) 
1 -7 

rexp(— jt) — 7exp(— Tt) 



1 



r- 7 



(El) 
(E2) 
(E3) 



(E4) 
(E5) 

(E6) 



The dynamics of a monomer is thus independent of the 
value of a. 

One also easily sees that the EST 77(a) is given by: 



77(a) 



1 
f 



1 

1 



(E7) 



Therefore, the transport efficiency always increases with 
increasing 7 and T. For example a high input with high 
output leads to a high efficiency and vice-versa. 
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